arXiv:1508.00398v3 [cond-mat.quant-gas] 3 Apr 2016 


Anomalous pseudogap in population imbalanced Fermi superfluids 


Madhupama Karmakar and Pinaki Majumdar 
Harish-Chandra Research Institute, Chhatnag Road, Jhunsi, Allahabad 211 019, India 

(Dated: April 5, 2016) 

In a Fermi superfluid increasing population imbalance leads initially to reduction of the transition temperature, 
then the appearance of modulated Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states, and finally the suppression 
of pairing itself. For interaction strength such that the ‘balanced’ system has a normal state pseudogap, increas¬ 
ing imbalance reveals anomalous spectral behavior. At a fixed weak imbalance (small polarization) the stable 
homogeneous superfluid occurs only above a certain temperature. The density of states has a minimum at the 
Fermi level, then a weak peak within the gap , and then the large, gap edge, coherence features. On heating, 
this non monotonic energy dependence changes to a more conventional fluctuation driven pseudogap, with a 
monotonic energy dependence. At large imbalance the ground state is FFLO and ‘pseudogapped’ due to the 
modulated order. It changes to a gapless normal state on heating, and then shows a pseudogap again at a higher 
temperature. These weak imbalance and strong imbalance features both involve effects well beyond mean field 
theory. We establish them by using a Monte Carlo technique on large lattices, motivate the results in terms of 
the pairing field distribution, and compare them to spectroscopic results in the imbalanced unitary Fermi gas. 


I. INTRODUCTION 

Attractive interaction between fermions lead to pairing and 
a superfluid or superconducting ground state. If the interac¬ 
tion strength is large the pairing effects show up at a high 
temperature, T pair , via a pseudogap in the density of states, 
while the transition to a superfluid occurs at a lower temper¬ 
ature T c . Between T pair and T c the pseudogap deepens and 
for T < T c one expects a full gap in the density of states. The 
loss in low energy spectral weight with reducing temperature 
is monotonic 1 . 

Since the pairing occurs between time reversed states k f 
and — k | an imbalance in the population of the ‘up’ and 
‘down’ species, with the total fixed, tends to reduce the pair¬ 
ing amplitude and the T c . The imbalance can be achieved by 
applying a magnetic field, as in the solid statJ 2 , or by loading 
different numbers of up and down fermions into a trap - as in 
cold atom experiments 3 5 . 

The effect of increasing field, h, or growing population im¬ 
balance, P, is qualitatively similar: it suppresses the T c of 
the homogeneous superfluid (SF), promotes a ‘phase modu¬ 
lated’ Fulde-Ferrell (FF) or an ‘amplitude modulated’ Larkin- 
Ovchinnikov (LO) state in the intermediate P regime, and fi¬ 
nally destroys pairing altogethei®. There are also crucial dif¬ 
ferences. For example, while there is a stable ground state for 
all applied fields , there is no stable ground state for polariza¬ 
tion 0 < P < P c i, where P c \ is the threshold polarization of 
the FFLO state. At zero temperature the system can either be 
an unpolarised (homogeneous) superfluid (USF) with P = 0, 
or a modulated state with P cl < P < P c2 , or a homogeneous 
‘normal’ state with P > P c2 . There is an unstable (phase sep¬ 
aration) window 0 < P < P c \. This window shrinks with 
increasing temperature vanishing at a tricritical point. The 
general features are well understood theoretically^ 121 and have 
been verified in cold Fermi gases at unitarity 313 . 

What is less well known is the changing spectral character 
of the imbalanced SF with increasing temperature. We estab¬ 
lish this for the two dimensional attractive Hubbard model at 
intermediate coupling, where the T c is highest, using a Monte 
Carlo approach. We discover the following: 


1. At low imbalance the homogeneous superfluid is stable 
only above some temperature T un (P) and for T > T un 
has a pseudogap with a strange non monotonic low en¬ 
ergy density of states. Heating this leads to a crossover 
to a more conventional pseudogap phase with a mono¬ 
tonic energy dependence. 

2. For P > P c i, where the ground state is FFLO, increas¬ 
ing temperature generates multiple crossovers: from the 
‘pseudogapped’ FFLO to a gapless normal state, and, 
beyond a higher temperature, to another pseudogapped 
normal state. 

3. The non paired ground state, with P > P c2 , is gap¬ 
less but on heating develops a weak pseudogap beyond 
a certain temperature. This phase is observed in a sit¬ 
uation where there is no order at any temperature, i.e, 
T c = 0. 


II. MODEL AND METHOD 

We study the attractive two dimensional Hubbard model on 
a square lattice in the presence of a magnetic field: 

H = H 0 - h^2a iz - \U\ y^n it n4 (1) 

i i 

with, H 0 = ~ ^u) c !a c P’ where Uj = -t only 

for nearest neighbor hopping and is zero otherwise. a iz = 
(1/2) (n^ — riii). We will set t = 1 as the reference energy 
scale, pt is the chemical potential and h is the applied mag¬ 
netic field in the z direction. U > 0 is the strength of on-site 
attraction. We will use U/t = 4 and p = —0.2 1 ( n ~ 0.94). 

We have discussed our method of solving the model in de¬ 
tail elsewhere 1415 so we touch on it only briefly. We use a 
‘single channel’ Hubbard-Stratonovich (HS) decomposition 
of the interaction in terms of an auxiliary complex scalar field 
A i(r) = \Ai(r)\e ldi ^ T \ This maps the original interacting 
problem to that of non interacting fermions in a space-time 



fluctuating field A*(r), the price to pay is the additional ‘av¬ 
eraging’ over all configurations of A*(r). Quantum Monte 
Carlo performs this averaging without further approximation, 
mean field theory (MFT) restricts A i(r) to a time independent 
spatially periodic function, and dynamical mean field theory 
(DMFT) treats A*(r) as a ‘single site’ time dependent func¬ 
tion A (r). 

We drop the ‘time’ dependence of A, but retain the full spa¬ 
tial dependence: A^(r) —» A^. This corresponds to dropping 
quantum fluctuations but retaining all the classical thermal 
fluctuations. Equivalently one can think of Fourier transform¬ 
ing the A i(r) into bosonic Matsubara modes A^(D n ) and our 
approach corresponds to retaining the £l n = 0 mode. This 
has been called a ‘static path approximation’ (SPA ) 16 17 to the 
functional integral for the partition function. 

SPA retains classical fluctuations of arbitrary magnitude, 
around the saddle point, but no quantum (fl n ^ 0 ) fluctua¬ 
tion. At T = 0, since the classical fluctuations die off, SPA 
reduces to standard Bogoliubov-de Gennes (BdG) mean field 
theory (MFT). For T 7 ^ 0, however, it considers not only 
the saddle point configuration but all configurations involving 
classical amplitude and phase fluctuations of the order param¬ 
eter. The BdG equations are solved in all these configurations 
to compute the thermally averaged properties. This approach 
suppresses the order parameter much quicker than MFT. Also, 
since the = 0 mode dominates the exact partition function, 
SPA becomes exact as T — 00. The approximation leads to 
the coupled equations: 

H e ff = Hq h ^ ^ (Ji z T~ ^ ^ (A^c^c]^ H- h.c ) H- H c i 

i i 

P{Ai} ex Tr C( 4 er ,3H ' !ff (2) 

where H c i = ^ is the stiffness cost associated with 

the now classical auxiliary field. The upper equation de¬ 
scribes fermions in a pairing background A i, while the lower 
equation defines the probability distribution P{A^} associ¬ 
ated with a {A^} configuration. 

We generate samples of { A^} using a Metropolis algorithm, 
diagonalising H e f f for the update cost. In order to make the 
study numerically less expensive the Monte Carlo is imple¬ 
mented using a cluster approximation 18 . After equilibriation 
at a given h, T we calculate the following: 

P(h,T) = ((l/N)J2(nit-n ii )) 

S< l (h,T) = ±J2( A i A*)e^-^ 

i,3 

N t (u) = {(l/N^Kl 2 ^ ~ E n )) 

i,n 

where P(h,T) is the polarization (not be confused with the 
probability distribution of the pairing field discussed later), 
S q (/i, T) is the pairing structure factor, and N^(uj) is the up 
spin density of states (DOS). The down spin DOS is symmet¬ 
rically shifted. u\ are the usual BdG eigenfunctions and E n 
are the BdG eigenvalues in an equilibrium configuration. An¬ 
gular brackets indicate thermal average. Notice that the DOS 



FIG. 1. The polarization (P) - temperature (T) phase diagram of 
the intermediate coupling attractive Hubbard model (U — At) and 
n ~ 0.94 in two dimensions. The thermodynamic phases include 
a homogeneous superfluid, FFLO states, a partially polarized Fermi 
liquid (PPFL), and an ‘unstable’ phase separated region. The spin 
resolved density of states can either have a hard gap, pseudogaps 
of different character, PG-a or PG-b (see text), or be gapless. The 
superfluid to PPFL transition is indicated as T c2 and T c i, for 2nd 
and 1 st order transitions respectively, and the unstable to superfluid 
boundary as T un . The PG-a to PG-b crossover within the superfluid 
is denoted as Tpgl and the gapless to pseudogap crossover at large P 
denoted as Tpg 2- The variation of T at fixed P can take the system 
through different pseudogap regimes. 


is also the k-summed spectral function A a ( k, uj). Formally, 


A a (k,uj) = — (1/ 7r)/mG cr (k, lc) 


with, G a (k,cj) = Um^o G a (k, iu3 n )\%uj n ^u)+i T ) 

where G a (k, iuj n ) is the imaginary frequency transform of 
(ckcr(^)cj ccr (0)). The BdG quasiparticles are related to the 
original fermions via the eigenfunction matrix so the Green’s 
function and spectral functions can be readily computed 19 . 

We wish to probe results at fixed P, for varying T, as 
would be the case in population imbalanced cold Fermi gases. 
Within the Hamiltonian formulation P is a derived quantity, 
dependent on ft, T and the other model parameters. We there¬ 
fore solve the h — T problem first and then construct a P — T 
phase diagram out of it. 

Since mean field theory is widely used to study the imbal¬ 
anced Fermi gas it is useful to point out that, apart from the 
gross overestimate of T c , mean field theory predicts that the 
homogeneous SF is always gapped , and the normal state al¬ 
ways gapless. The notion of a pseudogap does not figure in 

ifD. 








FIG. 2. Temperature dependence of the up spin DOS at (a) P ~ 0.010, (d) P ~ 0.075, (e) P ~ 0.300 and (f) P ~ 0.375. The plots are 
vertically shifted for clarity and the uj axis is shifted by +h to center the gap feature on the origin. In (a) the DOS starts with PG-a character 
(see text) at the lowest accessible temperature, transforms to PG-b on heating (all within the SF phase), then to a pseudogapped normal state. In 
(d) the low T pseudogap weakens but persists across the SF to PPFL transition and on to high T. In (e) the ground state is ‘pseudogapped’ due 
to FFLO order. This is lost on transition to the PPFL but a weak PG reappears at high T. In (f) the PPFL ground state is gapless but develops 
a PG on heating. Panel (b) shows the DOS behavior at P ~ 0.010, focusing on the low energy window. Notice the change in character with 
increasing T. Panel (c) shows the energy derivative of the DOS at P — 0.01, the multiple changes in sign at low T indicate the multiple 
minima-maximum structure. The sign changes occur only for the two lowest T, which we have called the PG-a phase. 


III. RESULTS 


1. Ground state: 


A. Phase diagram 


Fig.l shows the P — T phase diagram inferred from the 
h — T phase diagram (shown later, and also in ref. 15 ). While 
the thermodynamic phase boundaries are determined based on 
S'q(T), the various crossovers are determined based on the 
behavior of the DOS. A suppression in the density of states 
around the Fermi level (but not a hard gap) is characterized as 
a pseudogap (PG). In our case we classify the pseudogaps as 
PG-a and PG-b depending on the detailed spectral behavior in 
the subgap region. Considering the cc > 0 window, PG-a has 
a minimum at cc = 0, a local maximum at a low energy cc maa: , 
a second minimum at an energy cc m ^ n , say, and then the gap 
edge coherence features. PG-b involves a minimum only at 
cc = 0. We first recapitulate the ground state then move to 
thermal properties. 


In terms of an applied magnetic field, the ground state is a 
homogeneous unpolarised SF for fields h < h c \ ~ 0.9 1. At 
h c i it makes a transition to a finite polarization FFLO state, 
with P = P c i. The FFLO regime is from h c i to h c 2 , with 
the polarization growing from P c 1 to P c 2 , and beyond h c 2 the 
ground state has no pairing and is a partially polarized Fermi 
liquid (PPFL). 


In terms of polarization, at T = 0 the entire 0 < h < h c \ 
region, with P = 0, collapses to the origin. The polarization 
regime 0 < P < P c \ is unstable, with no homogeneous phase 
allowed, followed by the FFLO and PPFL at higher P. The 
P = 0 superfluid is gapped, the FFLO has ‘pseudogaps’, i.e, 
depression in its DOS due to the spatial modulation, while the 
PPFL is gapless. 




























2. Thermal properties: 


At finite T the homogeneous SF does allow finite polariza¬ 
tion and is called the ‘breached pair’ (BP) state 15 . It occupies 
an increasing P window, as T increases, at the expense of the 
unstable region. The temperature above which the SF is stable 
is T un (P). At an even higher temperature, T c2 , superfluidity 
is lost through a Berezinskii-Kosterlitz-Thouless (BKT) tran¬ 
sition. We infer the T c scale from So(h, T ). 

The upper right edge of the unstable region defines the 
scale, T c i , for a first order transition between the SF and the 
normal state. The FFLO has a very low T c and transits to a 
gapless normal state on heating. This gapless phase defines a 
large part of the high P low temperature window. 

The SF phase can be gapped or pseudogapped as Fig.l 
shows, the FFLO is pseudogapped, while the PPFL is gap¬ 
less at low T and pseudogapped at high T. We will discuss 
the origin of these behavior in the discussion section, but the 
phase diagram indicates that varying temperature at a fixed 
polarization can lead to multiple changes in spectral charac¬ 
ter. This is directly relevant for cold Fermi gases where one 
works at fixed imbalance rather than a fixed applied field®. 


B. Density of states 

Fig.2 shows the spin resolved DOS for fixed P cross sec¬ 
tions through the P — T phase diagram in Fig.l. Two of the 
panel, (a) and (d), traverse the BP part of the phase diagram, 
(e) starts with a FFLO ground state, while (f) has a PPFL 
ground state. The ground states vary widely, we discuss the 
cases one by one. 

Fig.2(a) is for P = 0.01. The T = 0 state in this case would 
be in the unstable window and stable phases are defined in this 
case only for T > 0.035L All the way from this temperature 
to T c rsj 0.1 It the system is a BP superfluid. The behavior 
of the DOS, however, changes multiple times. For 0.035T < 
T < 0.08£ the (shifted) DOS has its absolute minimum at 
uj = 0 but then a local maximum at a small scale c Umax (see 
plot) and a local minimum after that at uj min and then the 
sharp rise at the gap edge. While subgap density of states 
by itself is not surprising this additional feature is unusual. 
Increasing T beyond ~ 0.08T leads to a more common PG 
with a monotonic increase in TV (a;) as cc increases. 

Fig.2(b) focuses on the low energy behavior of the data in 
2(a). The lower scales clearly reveal the non monotonic DOS 
at low T. We fitted this behavior to an approximate form 
N(uj) = a + buo 2 + ccc 4 + dec 6 , and, after extracting the coef¬ 
ficients a, 6, c, d, plotted the derivative dN/duo in Fig.2(c). 
The multiple zero crossings in the low T data highlight the 
multiple extrema. 

The result in Fig.2(d) is at P = 0.075. The stable BP win¬ 
dow at this P starts at T ~ 0.07 1. Heating this BP state leads 
to a 2nd order transition at T c ~ 0.14L All the way from 
T ~ 0.07 1 through T c to high temperature in the PPFL phase 
the system is pseudogapped, with the character that we called 
PG-b, with monotonic N(u). The coherence peaks vanish, 



FIG. 3. Magnetic field-temperature (h — T ) phase diagram showing 
the varying h(T) that is needed to maintain the polarization fixed 
at the values chosen in Fig.2. In general when pairing effects are 
strong and the low T phase is a BP superfluid the magnetic field has 
to sharply rise with reducing T to maintain P fixed. This is clearly 
visible for P — 0.01 and weakly so for P — 0.075. In the higher 
field regime, where the low T state is FFLO or PPFL the magnetic 
field needs to increase somewhat with increasing T, as in a normal 
metal, to maintain a fixed P. The shaded region within the BP phase 
roughly correspond to the anomalous PG phase (see text for details). 


the depth of the pseudogap lessens, and the low energy weight 
increases with increasing T. 

Fig.2(e) is for P ~ 0.30. The ground state in this case is an 
axial stripe LO phase, with Q = (7r/3, 0). The periodic mod¬ 
ulation of A i in the low T state leads to multiple bands and a 
redistribution of the tight binding spectral weight. The pres¬ 
ence of a finite Q modulation leads to pairing of fermions be¬ 
tween |kt) and | — k±Q|) states (rather than the time reversed 
|kt) and |—k|)) and give rise to multiple branches in the dis¬ 
persion. Some of these branches cross the Fermi level leading 
to a finite but depressed DOS around uj ~ 0. This is not a 
fluctuation induced pseudogap but a ‘band structure’ effect. 
The T window immediately above T c , T ~ [O.Olt, 0.05t], 
in Fig.2(e), is gapless. The pairing field is small and disor¬ 
dered and does not affect the DOS noticeably. However, at the 
two highest temperatures in Fig.2.(e) a weak pseudogap again 
emerges. This is due to the thermally induced growth in the 
mean magnitude (|A*|) with the phase variables 0i remaining 
random. At this polarization one has a pseudogap-gapless- 
pseudogap re-entrance with increasing temperature. 

Finally, Fig.2(f) at P = 0.375 has a non ordered ground 
state with no \ at any site. This is a system at strong in¬ 
teraction, U = 4 1, where the imbalance suppresses pairing at 
T = 0. As a result the low temperature phase is gapless. With 
increasing T the (| A^ |) grows quickly, as in the FFLO win¬ 
dow, and leads to a pseudogap for T > 0.07L This is a simple 
instance of a pseudogap emerging without the presence of any 
order in the low temperature state. 
















FIG. 4. The T dependence of the pairing field distribution for different P. The x axis is normalized with respect to the T = 0 unpolarised SF 
amplitude. The panels are for ( a) P ~ 0.01, (b) P ~ 0.075, (c) P ~ 0.300 and (d) P ~ 0.375. 


IV. DISCUSSION 

We first discuss the origin of the effects seen in Fig.l and 
Fig.2, and then earlier theory efforts and experimental data on 
this topic. 


A. h — T trajectory 

Fig.3 shows the field-temperature phase diagram from 
which the P — T phase diagram is derived. The basic phases 
have been discussed earlier in ref. 15 . It is included here mainly 
to show the rather unusual h — T paths that lie behind the con¬ 
stant P cross sections in Fig.l. 

At P = 0.01 the low T part of the trajectory passes close to 
the BP-LO phase boundary. We explored the h — T neighbor¬ 
hood and discovered that in the window shaded grey the DOS 
indeed displays the peculiar features observed in Fig.2(a). The 
proximity to the BP-LO phase boundary, at finite T, suggests 
that while much of the A* would have (large) values appro¬ 
priate to the BP phase, there could be a small fraction of sites 
that have the lower values appropriate to the LO phase. We 
will examine this aspect in the distributions later. 


The P ~ 0.075 trajectory transits from a high T BP phase 
to the PPFL. The slope of the h(T) curve is much smaller than 
for P = 0.01 due to the weaker suppression of Pauli suscep¬ 
tibility. The P ~ 0.30 system passes from a low T LO phase 
to a PPFL, while the P ~ 0.375 path is always in a unordered 
phase but displays a high T pseudogap due to thermal fluctu¬ 
ations. For these two cases, where the ground state does not 
have Pauli susceptibility suppression, the h(T) curves gener¬ 
ally have a weak positive slope. In this window the ‘myste¬ 
rious’ effect is the presence of a high T pseudogap above a 
gapless low temperature state. Addressing this requires a look 
at the pairing field distributions. 


B. Pairing field distribution 

Fig.4 shows the distribution of the pairing field for the dif¬ 
ferent fixed P cross sections shown in Fig.2. Fig.4.(a) shows 
the pairing field distributions for DOS shown in Fig.2.(a), and 
so on for the other panels. 

The anomalous low energy features in Fig.2(a), we believe, 
arise from the relatively large number of small A i sites in a 
typical configuration in this part of parameter space. This 
arises due to the proximity to the low A FFLO state. While the 






























FIG. 5. Size dependence of the DOS at a polarization P = 0.01 and three temperatures. The curves are y shifted by equal amounts for clarity. 
We observe that (i) L — 24 and L — 32 results are quite similar (and different from L — 16), and (ii) the DOS at low temperature, panels (a) 
and (b), have a consistent non monotonic low energy feature. 


mean A is ~ A 0 (the T — 0 mean field value), at T ~ 0.04£ 
in Fig.4(a), if the system were fully in the LO regime, Fig.4(c), 
the mean A would be ~ 0.5Ao, see the plot for T = 0.05L 
We think that proximity to the LO phase boundary makes the 
system have a small fraction of small A sites whose magni¬ 
tude is ^ 0.5 the typical value. As a result the DOS generates 
weight at an energy that is roughly 0.5 the gap edge 

value. Since the P = 0.01 trajectory is the only one that bor¬ 
ders the BP-LO boundary the effect is not visible at other P. 

Fig.4.(b) shows a monotonic increase in the width of 
P(|A^|) with increasing T, and the increasing ‘disorder’ in 
the pairing field steadily weakens the PG in the DOS and in¬ 
creases the low energy spectral weight. In Fig.4.(c) the lowest 
T shows the ‘multimode’ distribution of the amplitude mod¬ 
ulated LO state. At this, and the next higher T, the LO mod¬ 
ulations lead to a ‘pseudogapped’ phase but with the loss of 
order at the first order transition, T ~ 0.005£, the mean am¬ 
plitude collapses. A gapless state emerges, and survives to 
some larger temperature, where the thermal fluctuations have 
created enough large A^ sites to generate a weak pseudogap. 
The behavior in 4.(d) is essentially similar to 4.(c), with the 
LO part removed: in this case one just has a gapless to weak 
PG crossover. 


C. Earlier work 

Much of the earlier work on pseudogaps in imbalanced su¬ 
perfluids is in the continuum context. A thermodynamic re¬ 
entrance has been pointed out in the continuum unitary gas 9 10 
at low P. The authors there had carried out an extensive in¬ 
vestigation in terms of interaction strength (l/kpa), polar¬ 
ization, and temperature and observed that near l/kpa = 0 
the low P s uperfl uid state has an upper and lower bound in 
temperatur j 9 l 1Q l 2Q l At T = 0 there is no stable superfluid, 
while an intermediate SF phase is realized at T / 0. The su¬ 
perfluid is gapped (by construction in this theory) the higher 
and lower T ‘normal’ states are pseudogapped. The present 
theory allows for pairing field fluctuations, consequently we 
find that within the SF phase the spectrum is not always 


gapped out completely and there are finite low energy weight 
even at the lowest accessible temperature. Despite the differ¬ 
ent model, a quite different approximation, and a thermody¬ 
namic stability (rather than spectral change) argument offered 
by the authors, we think the similarity with our result is not 
coincidental. The changing spectral character at the small po¬ 
larization regime bears closer experimental search. 

The spectrum in imbalanced superconductors can be 
probed by tunneling but we are not aware of such stud¬ 
ies. Spectral studies are more visible in cold atomic super¬ 
fluids, although the results are complicated by the effect of 
a trap. The experiments measure a ‘radio frequency’ (RF) 
current 21 2 3 rather than the density of states. However, the 
RF response I a (uj) is a weighted sum 2425 of the same spec¬ 
tral function A a (k, ut) = — (l/7r)/mG cr (k, uj) that defines 
our density of states N a (cj). 

Experiments on the unitary Fermi gas indicate that the T c 
falls from ~ 0.2 T F at P = 0 to zero at P c ~ 0.75. The 
‘pair formation’ scale T pa i r , however, seems to be 0(Tp ) 
even when P crosses P c and T c falls to zero 13 . These exper¬ 
iments suggest that a PG is obtained over a wide temperature 
window all the way from P = 0 to the highest P. Our results 
agree with this, except at the high P end where we observe a 
gapless low temperature state. This difference arises from our 
neglect of quantum fluctuations in the A i(r), which can gen¬ 
erate a non trivial Fermi liquid ground state 2 ^U, preserving 
pairing but suppressing condensation. 


D. Finite size effect 

The results discussed in this paper correspond to a lattice 
size of L = 32. The system size is reasonable in comparison 
to the existing literature in which the thermal physics of spin 
imbalanced system are being explored^. In order to verify the 
robustness of the spectral features, we have further computed 
our results on various other system sizes. Fig.(5) shows the 
spin resolved DOS calculated at three different lattice sizes of 
L = 16, 24 and 32, for a fixed polarization, highlighting three 
different temperature regimes. We observed that the anoma- 






























lous behavior of the spin resolved DOS is fairly robust and 
persists even at L = 24, and is not an artifact of finite size 
effect. 


V. CONCLUSION 

We have established the thermal phase diagram of popula¬ 
tion imbalanced lattice fermions near the BCS-BEC crossover 
and observed that although the thermodynamic phases exhibit 
expected behavior with temperature the spectral features are 
very counter-intuitive. The low temperature superfluid at a 
fixed small polarization has a pseudogap with a distinct max¬ 


imum in the density of states in the subgap region, crossing 
over to a more conventional featureless pseudogap with in¬ 
creasing temperature. The unusual low (but finite) tempera¬ 
ture result arises due to LO like fluctuations in the BP phase 
near the BP-LO boundary. At large polarization, in the LO 
phase, the system undergoes a pseudogap to gapless to pseu¬ 
dogap crossover with increasing temperature, At even larger 
polarization, where the ground state is gapless and has no pair¬ 
ing and long range order, heating the system leads to pseudo¬ 
gap formation above a certain temperature. All these effects 
are a consequence of thermal fluctuations beyond the mean 
field scheme usually used to analyze these models. 

We acknowledge use of the HPC clusters at HRI and thank 
Nyayabanta Swain for comments. 
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